Be(p,n)i^B reaction at 69 MeV in inverse kinematics 
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A Gamow- Teller (GT) transition Irom the drip-line nucleus ^'*Be to '^^B was studied via the (p,n) 
reaction in inverse kinematics using a secondary ^'*Be beam at 69 MeV/nucleon. The invariant 
mass method is employed to reconstruct the energy spectrum. A peak is observed at an excitation 
energy of 1.27(2) MeV in "B, together with bumps at 2.08 and 4.06(5) MeV. The observed forward 
peaking of the state at 1.27 MeV and a good description for the differential cross section, obtained 
with a DWBA calculation provide support for the l""" assignment to this state. By extrapolating 
the cross section to zero momentum transfer the GT-transition strength is deduced. The value is 
found to compare well with that reported in a /3-delayed neutron emission study. 

Keywords: Gamow- Teller transition, ^H(^'*Bc,'^^B-(-n), Unbound states. Distorted-wave approximation 



The Gamow- Teller (GT) transition, characterized by 
an angular momentum transfer of AL=0, a spin transfer 
of AS'=1, and an isospin transfer of AT=1 between ini- 
tial and final nuclei, represents one of the fundamental 
modes of nuclear excitations. Its properties are impor- 
tant not only for understanding spin-isospin correlations 
in nuclei [1], but also for their implications for a variety 
of astrophysical phenomena where weak processes play 
a key role, such as late stages of stellar evolution and 
neutrino nucleosynthesis [2]. The rapid increase in the 
availability of Radioactive Isotope (RI) beams in recent 
years has made it possible to study GT strengths in nuclei 
far from the /3-stability line. Although /3 decay has been 
the major source of information about such strengths, 
it might be expected that nuclear charge-exchange (CE) 
reactions play an important role [3], with their distin- 
guishing feature of being unconstrained by any lim- 
itations. We performed an experimental study of the 
(p, n) reaction on a neutron-rich nucleus ^"^Be using the 
invariant mass method in inverse kinematics. The data 
on the ^^Be(g.s.) — >-^^B(l^) transition process enabled 
the extraction of i3(GT) from forward angle (p, n) cross 
sections. Studies were also made on states populated in 
the iH("Be,"B+n) reaction. 

The -'^■^Be nucleus is located at the neutron drip line. 
A two-neutron halo structure has been suggested by in- 
clusive [4, 5] as well as exclusive [6] reaction studies. Un- 
ambiguous identification of the first 2+ state was made 
recently [7]. The /? decay of ^^Be is known to be dom- 
inated by one-neutron emission with a branching ratio 
compatible with 100% [8]. A concentration of the decay 



strength in a state at 1.28 MeV, which is located just 
above the one-neutron separation energy of S'„=0.97(2) 
MeV [11] in ^'^B, was found in energy spectra of decay 
neutrons [9, 10]. The smah log ft value of 3.68(5) [10] 
has indicated that the decay is due to a GT transition, 
providing the 1^ assignment to this state. 

Existing CE studies on /3-unstable nuclei have been 
based on two major methods: the first involves final 
states in a bound region and the missing mass technique 
was applied, where momenta of heavy charged ejectile 
were measured either by using a spectrograph in the 
6He(p,n)6Li [12, 13] and 34p(7Li/Be-h 7)^*81 [14] reac- 
tions, or a telescope in the ^He(p, n)^Li reaction [15]. 
The second method involves unbound final states, and 
the invariant mass technique was adopted. All particles 
from a decay of an unbound state were detected using a 
plastic counter hodoscope in the ^^Li(p, n)^^Be [16] and 
^''Be(p, n)^^B [17] reactions; isobaric analog states in the 
residual nuclei were identified. Until now, there exist 
only a limited number of CE studies involving unstable 
nuclear beams; it is desirable to further explore experi- 
mental approaches which permit reliable determination 
of relevant nuclear structure information. The present 
study, succeeding to the latter works, represents the first 
extraction of the GT strength on an unstable nucleus, 
leading to a particle unbound state, in an inverse kine- 
matics CE reaction. Independent detectors were used for 
different decay products, charged fragment and neutron, 
to achieve high enough acceptance near the threshold, 
together with a liquid hydrogen target to increase the re- 
action yield. In inverse kinematics, decay neutrons emit- 
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ted from the unbound system, which is Lorentz boosted 
with a beam-hke velocity in the laboratory, acquire high 
energies and can be detected with high efficiencies. With 
this new feature, it is shown that the inverse kinematics 
CE reaction accompanied by neutron decay provides an 
alternative approach for nuclear structure, to /3-delayed 
neutron spectroscopy. 

The extraction of the GT strength in CE reactions re- 
lies on the proportionality relationship between the re- 
duced transition probability B{GT) and the cross section 
at zero momentum transfer (q=0) [18, 19]: 
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ND\Jar\^B{GT). (1) 



Here is the reduced energy in the initial (final) 

channel and {J^tI is the volume integral of the cen- 
tral part of the effective nucleon-nucleon {NN) interac- 
tion in the spin-isospin channel. The distortion factor 
ND = {da/dn)^^o/{da/dn)^^Q is the ratio of distorted- 
wave (DW) and plane-wave (PW) cross sections. The 
cross section {da / dVL)^^^ , defining the L=Q contribution 
at g=0, is obtained by extrapolating the 0° cross sec- 
tion to g=0 using the ratio of distorted-wave Born ap- 
proximation (DWB A) cross sections at g=0 and 0° , with 
an additional factor given by the ratio of the cross sec- 
tion calculated with only the central terms of the NN 
interaction and that with the full force components. Al- 
though the above relationship [Eq.(l)] has been mainly 
used for CE reaction data at bombarding energies above 
100 MeV/nucleon, it is expected to be applicable for 
(p, n) scattering data involving transitions with appre- 
ciable strengths, at energies of the present work (around 
Ep=7Q MeV), due to the simple reaction mechanism of 
the nucleon scattering and to the availability of modern 
nucleon- nucleus optical model potentials (OMPs). 

The experiment was performed at the projectile- 
fragment separator RIPS [20], which is a part of the 
RIBF accelerator complex operated by RIKEN Nishina 
Center. The setup was identical to the one used in 
Ref. [21]. The radioactive beam of ^''Be was produced 
from a 100-MeV/nucleon ^^O beam, which impinged on 
a 6-mm-thick Be target. The typical beam intensity was 
7 kcps with a momentum spread of AP/ P—zL2%. The 
beam profile was monitored by a pair of parallel-plate 
avalanche counters (PPACs) , placed upstream of the tar- 
get. The secondary target was a liquid hydrogen tar- 
get [22] with a cylindrical shape: 3 cm in diameter and 
229±6 mg/cm^ in thickness. The average energy of ^^Be 
at the target center was 69 MeV/nucleon. The target 
was surrounded by a Nal(Tl) scintillator array to detect 
7 rays from charged fragments. The fragments bent by 
a dipole magnet behind the target were detected by a 
plastic counter hodoscope, which gave Z identification 
by AE. Their trajectory was measured by a set of multi- 
wire drift chambers placed before and after the magnet. 
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FIG. 1. Relative energy spectra of the ^H(^''Be,"B + n) re- 
action for (a) 6lcm=0°-16° and (b) 6lcm=16°-48°. 



which, together with time of flight, gave mass identifica- 
tion. Neutrons were detected by two walls of plastic scin- 
tillator arrays placed 4.6 and 5.8 m downstream of the 
target. The neutron detection efficiency was 24.1 ± 0.8% 
for a threshold setting of 4 MeVee. The relative energy 
Eic\ of the final system was calculated from the momen- 
tum vectors of the charged fragment and the neutron. 
The bending magnet for the fragment allowed the two 
detector systems to be located separately. This, together 
with the kinematical focusing of the reaction products 
toward forward angles, ensured that the phase space vol- 
ume of the final decaying system was covered efficiently. 

Figure 1 shows relative energy spectra obtained at two 
different angular bins: (a) 0cm=O°-16° and (b) 0cm=16°- 
48°. The effects of finite detector acceptances were cor- 
rected. Background contributions due to window mate- 
rials, measured with an empty target, were less than 0.1 
mb/MeV, and were subtracted in the figure. Error bars 
are statistical ones. A sharp resonance peak, evident at 
-Erei^O.S MeV, is most prominent at forward angles. A 
bump structure at iJidWl.l MeV is visible only at back- 
ward angles, and might correspond to the reported 4~, 
2.08 MeV state observed in the "c(7Li,7Be)"B [23] and 
i2C("C,i2n)"B [24] reactions. Since this state is weakly 
populated, the energy of 2.08 MeV was adopted in the 
subsequent analysis. Another bump at E,-ci~S MeV was 
observed for the first time in this study. For the three 
states no coincident 7 rays were identified. 

The spectra were analyzed using single Breit-Wigner 
shape functions as described in Ref. [25] to extract the 
resonance energy Er, width F,., and peak yield. The line 
shape was obtained in a Monte Carlo simulation which 
includes the detector response, beam profile, and target 
thickness. The relative energy resolution was well de- 
scribed by the relationship: A 25^1=0.20^/^x^1 MeV(rms). 
For strengths beneath the resonance, thermal emission of 
a neutron from the continuum [26] was assumed. A con- 
stant background, which simulates various effects such 
as overlapping broad resonances and decay to an excited 
residue was also introduced when needed. Solid lines in 
Fig. 1 are the results of the fit; the assumed continuum- 
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TABLE I. Resonance parameters in the present analysis. The 
excitation energy of Ex—2.08 MeV is from Refs. [23, 24]. 



Er (MeV) 


(MeV) 


Er (MeV) 


I (h) J" 


0.304(4) 


1.27(2) 


0.16(2) 


1 1+ 


1.11 


2.08 




2 4" 


3.09(5) 


4.06(5) 


[1.0(3), 1.2(5)] 


(1,2) (3+,3-) 



like strengths are indicated as dashed hnes. Resulting 
resonance parameters are summarized in Table I. The 
orbital angular momenta (l) of decay neutrons are also 
shown. The excitation energy {Ex=Er+Sn) of 1.27(2) 
MeV for the peak at i?rci~0.3 MeV agrees well with the 
reported value of 1.28(2) MeV for the first 1+ state in 
[9, 10]. The width of rr=0.16(2) MeV is consistent 
with the single-particle estimate of Fsp =0.31 MeV for 
an 1=1 resonance calculated in a Woods-Saxon potential, 
while it is significantly larger than r'sp=0.02 MeV for a 
d-wave, supporting the assumption of a p-wave decay. It 
is to be noted that the single particle width, Egp, is an 
estimate of the total decay width derived within the one- 
particle model in a simple nuclear potential, and that it 
provides an upper limit for the actual nucleus case where 
the configuration mixing is generally important. For the 
2.08 MeV state, the width was not reliably extracted due 
to limited signal-to-background ratio near the resonance; 
it was assumed to be dominated by the resolution. Note 
that an experimental upper limit of rr=0.l MeV [24] has 
been reported. For the 4.06 MeV state two different de- 
cay angular momenta, 1=1 and 2, were considered. The 
deduced widths are shown in Table I. The resonance en- 
ergy was similar in the two cases and varied within the 
errors. 

Figure 2 shows differential cross sections leading to the 
observed states. These were obtained by subdividing the 
acceptance solid angle to have a spectrum at each an- 
gular bin, and applying the fitting procedure to extract 
the peak yields. Error bars are statistical ones; system- 
atic uncertainties are estimated to be 6%, including those 
in the target thickness, neutron detection efficiency, and 
fitting procedure. The resolution in the scattering angle 
was estimated with the Monte Carlo simulation to vary 
from A6icm=2.3° in rms at 6'cm=4° to 4.5° at 6'cm=44°. 

Microscopic DWBA calculations were performed to 
clarify the nature of states observed and to have a means 
to extrapolate the 1"*" cross section to g=0. The OMP 
was taken from the semi-microscopic parametrization of 
JLMB [28]. Within the local density approximation, the 
local value of the finite nucleus OMP was deduced from 
the nuclear matter OMP. The proton and neutron radial 
density distributions were calculated in a mean field using 
the SkX interaction [29]. The effective NN interaction 
was the M3Y interaction [30] (force components labeled 
1,4,11,14,16,17, and 20 were used). The transition den- 
sity was obtained from the shell model using the code 
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FIG. 2. Angular distributions of the differential cross sec- 
tion leading to the three states observed (open circles). Data 
for the ^2C(p,n)^^N(l+,g.s.) reaction at 61.9 MeV [27] (open 
squares) are also shown. Normalization factors applied to the 
DWBA curves (solid lines) are indicated in the figure. Addi- 
tional curves, dashed and dotted lines, given for 0^— 7>1^ tran- 
sitions are DWBA results obtained, respectively, with only the 
central and tensor terms of the NN interaction. 



OXBASH [31]. The calculations used the WBT Hamil- 
tonian [32] in the spsdpf model space within the Ohuj 
configuration. The single-particle wave functions were 
generated in a harmonic oscillator well with an oscilla- 
tor parameter 6=2.04(6) fm, chosen to reproduce the rms 
matter radius of ^^Be [33]. The code DW81 [34] was used 
for the calculations. 

The DWBA cross sections are shown as solid lines in 
Fig. 2. Normalization factors required for the theoreti- 
cal curves to reproduce the data are also indicated. For 
the calculation of the 1.27 MeV state, the wave function 
of the first l"*" shell-model state predicted at 1.01 MeV 
was used, while for the 2.08 MeV state that of the first 
4~ state at 1.75 MeV was used. The measured cross 
sections leading to these states, which exhibit character- 
istic angular distribution patterns, are well reproduced 
by the calculations, corroborating the earlier spin-parity 
assignments of 1+ for the 1.27 MeV state [10] and 
A- for the 2.08 MeV state [23, 24]. For the 4.06 MeV 
state, the calculation was repeated by assuming differ- 
ent shell-model wave functions with various values. 
The angular distribution was found to be best described 
by the curve assuming the second 3"*" state predicted at 
4.20 MeV (red line). Another possible candidate was the 
second 3~ state at 5.70 MeV. Although the calculated 
cross section, shown in Fig. 2 as a black line, is slightly 
backward peaked, it is compatible with the data and the 
possibility could not be completely ruled out. With this 
a tentative J"^ assignment of 3"^ or 3~ is made to this 
state. It would be interesting if we could distinguish the 
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two cases by experiments with different sensitivity from 
the present one. 

The extraction of S(GT) leading to the 1+, 1.27 MeV 
state was done by using Eq.(l). The DWBA frame- 
work as described above was utiUzed to calculate the pa- 
rameter for momentum extrapolation and the distortion 
factor. In order to sec the contribution from the non- 
central terms of the NN interaction, two kinds of ad- 
ditional DWBA calculations were performed: one with 
only the central terms and the other with only the ten- 
sor terms, and the results of these calculations are shown 
as dashed and dotted lines in Fig. 2, respectively (the 
contribution from the spin-orbit terms turned out to be 
negligibly small). The ratio of the cross section calcu- 
lated with only the central terms to that with full force 
components was 0.961 at 9=0°, while it was closer to 
unity, 1.008, at q=0. Two additional sets of OMP pa- 
rameters [35, 36] were used to estimate the uncertainty 
due to the choice of the optical model. For calculations 
of the cross section at g=0, the distorting potential in 
the exit channel was re-evaluated at a neutron energy 
which satisfies the kinematical constraint. The volume 
integral iJar] of the NN interaction can be deduced effec- 
tively by analyzing the existing {p,n) 0+^1+ transition 
data at a similar incident energy, for which the B{GT) 
value is well established. The ^^C(p, n)^^N(l+, g.s.) re- 
action data at 61.9 MeV [27], with the corresponding 
B{GT) value of 0.875(6) [37], were used for this pur- 
pose. In the DWBA analysis, three sets of OMP param- 
eters [28, 35, 36] and the P(10-16)T shell-model wave 
function [32] were used. The oscillator parameter chosen 
was 6=1.87 fm [38], which reproduces the maxima of the 
electron scattering form factors. A result of the fit of the 
theoretical curve, calculated with the JLMB [28] OMP, 
to the data is shown as a solid line in Fig. 2. Dashed and 
dotted curves are again DWBA cross sections calculated 
with only the central and tensor terms of the A'^A'^ inter- 
action, respectively. For this transition, switching off the 
tensor terms had an effect to reduce the cross section by 
1.3% (-0.8%) at 0=0° ((7=0). By using Eq.(l), | J<,^| at 
61.9 MeV was deduced to be 212.1±3.1(stat)±9.3(syst) 
MeVfm^, where the systematic error includes the uncer- 
tainty due to the choice of optical parameters (2%). The 
value turns out to be consistent with the tabulated ones 
for the M3Y interaction [39]: ] J^t-|=211(247) MeVfm^ at 
£'p=61(40) MeV. From the proportionality between the 
cross section at q=0 and B{GT) [Eq.(l)], which assumes 
IJarl at 69 MeV obtained by making energy extrapola- 
tion using values for the M3Y interaction, the _B(GT) 
for the transition to the 1+ state at 1.27 MeV in ^''B 
was found to be 0.79±0.03(stat)±0.09(syst). The sys- 
tematic error includes the uncertainty due to the choice 
of optical parameters (4%). The value is in good agree- 
ment with the ^-decay value of B(GT)=0.80±0.09 [10], 
calculated from the ft value using the relationship: 
B{GT)=igA/gv)-'{6U7/ft) [1]. 



The shell-model prediction obtained using the WBT 
interaction [32] and the effective GT operators [40], 
within the Ohu basis, is B(GT)=0.88, which is consis- 
tent with the data. This suggests that the ground state of 
^^Be is predominantly of Ohoj nature. Precision data and 
improved knowledge on effective GT operators will allow 
further detailed analyses of the admixtures of higher-Tio; 
intruder components [41], in terms of B{GT). 

The ^^Be(g.s.) -^^^B{1^) transition exhausts only a 
minor fraction (-4%) of the GT sum rule: 3{N-Z)=18. 
Major GT strengths are expected to be located at previ- 
ously unexplored high excitation energy regions [42, 43]. 
CE studies probing such strengths arc a challenge (this 
might need either multiple-neutron detection within the 
current invariant mass scheme or recoiled-neutron detec- 
tion in the missing mass scheme) , but will shed new light 
on our understanding of spin-isospin correlations in nu- 
clei at extreme conditions of isospin and binding energy. 
The present study, quantifying the accuracy for a pro- 
cedure to extract i?(GT) from forward angle (p, n) cross 
sections on radioactive beams to be of order 15%, paves 
the way for future investigations in this direction. 

In summary, the {p, n) reaction on ^^Bc was measured 
at 69 MeV/nucleon using the invariant mass method. 
In decay energy spectra three resonance states were ob- 
served. A comparison of measured cross sections with 
DWBA calculations yielded a confirmation of the ear- 
lier assignments of l"*" and 4~ for the 1.27 and 2.08 
MeV states, respectively. It also led to a tentative J'^ 
assignment of 3^ or 3^ for a state newly observed at 
4.06 MeV. By extrapolating the cross section to q=0, the 
i?(GT) leading to the 1+ state was deduced, which com- 
pared well with the /?-decay value, providing the first 
successful demonstration of an experimental capability 
of extracting GT strengths on unstable nuclei via decay 
spectroscopy following CE reactions in inverse kinemat- 
ics. 
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